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ABSTRACT 

The  existence  of  travelling  wave  solutions  of  a  diffusion 
reaction-system  is  studied  via  the  generalized  Morse  index  of 
isolated  invariant  sets.  This  index  theory  is  analogous  to  degree 
theory,  and  the  method  of  proof  follows  lines  familiar  from  the 
latter  theory.  The  equations  in  question  are  "deformed"  to  a 
"standard"  system  where  the  index  can  be  easily  computed,  and  the 
existence  theorem  follows  from  the  "non-triviality"  of  the  index. 

The  index  theory  has  been  described  in  other  papers;  here  the 
main  job  is  to  construct  "isolating  neighborhoods"  which  are  analog¬ 
ous  (in  the  degree  theory)  to  open  sets  with  no  critical  points  on 
the  boundary.  Some  novel  means  of  locating  such  neighborhoods  are 
described. 

The  main  theorem  concerns  a  case  where  the  reaction  system  is 
a  competitive  system:  that  is,  the  growth  rate  of  one  population 
decreases  as  the  other  population  increases.  In  a  plausible  class 
of  such  models,  each  system  admits  three  stable  equilibria:  the 
two  at  which  one  population  is  at  a  maximal  stable  level  and  the 
other  is  eliminated,  and  one  at  which  both  populations  are  at  the 
zero  level.  In  general  there  are  (two)  travelling  waves  connecting 
one  or  the  other  of  the  first  two  equilibria  to  the  third.  The 
theorem  gives  a  criterion,  in  terms  of  the  relative  velocity  of 
these  two  waves,  in  order  that  there  be  a  third  wave  running  between 
the  first  two  equilibria. 
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SIGNIFICANCE  AND  EXPLANATION 


Two  populations  are  in  competition  if  an  increase  in  one  causes 
a  decrease  in  the  growth  rate  of  the  other.  In  the  absence  of  any 
other  effects,  the  interaction  of  the  two  populations  can  be  modeled 
by  a  system  of  ordinary  differential  equations  describing  the  change 
in  population  levels.  Typically  this  system  admits  three  attracting 
rest  points:  two  at  which  one  population  has  reached  its  maximum 
(stable)  level  and  the  other  has  been  eliminated,  and  the  third 
where  both  populations  are  at  zero  level. 

Now,  allowing  for  spatial  effects,  one  can  consider  a  (symmetric) 
situation  where  one  population  is  at  its  maximal  level  in  one  part 
of  the  space  (with  the  other  eliminated)  and  the  reverse  situation 
holds  in  the  other  part.  If  the  populations  are  able  to  move  in 
space,  it  can  be  expected  that  one  or  the  other  of  them  might  "win 
out"  over  the  whole  space.  This  report  concerns  mathematical  problems 
related  to  the  above  situation  where  the  spatial  movement  is  modeled 
by  diffusion. 

Among  a  class  of  mathematical  models  of  a  plausible  type,  two 
possibilities  are  indicated.  In  the  first,  one  of  the  two  pop¬ 
ulations  will  eventually  dominate  and  the  other  will  be  eliminated. 
This  corresponds  to  the  existence  of  a  "travelling  wave"  solution 
(of  the  diffusion-reaction  system)  which  connects  the  first  two  of 
the  above  mentioned  rest  points.  In  the  second  situation  this 
solution  does  not  exist;  instead  there  exists  a  "composition"  of 
two  travelling  waves,  each  of  which  connects  one  of  the  above  two 
rest  points  to  the  third  at  which  both  populations  are  zero. 

Evidently,  in  this  case  neither  population  is  strong  enough  to 
survive  the  conflict.  (However,  this  is  not  proved  in  this  report.) 

The  main  result  gives  a  simple  way  to  distinguish  these  two 
possibilities  on  the  basis  of  the  "relative  strengths"  of  the 
attractors. 


The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MRC,  and  not  with  the  authors  of 
this  report. 
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TRAVELLING  WAVE  SOLUTIONS  OF  A  COMPETITIVE  j  -  '  ■  >or 

REACTION-DIFFUSION  MODEL  ;  ' 


C.  Conley  and  R.  Gardner 


§  1.  Introduction 
A.  The  Equations 

This  report  concerns  the  ordinary  differential 
wave  solutions  of  a  diffusion-reaction  system.  The 
as  a  model  for  interspecific  competition  in  which  the 
over  a  "one-dimensional"  spatial  domain. 

Without  diffusion,  the  model  takes  the  form 


equations  for  travelling 
latter  can  be  interpreted 
two  species  diffuse 


(1) 


U1  =  ulFl(ul*u2> 

u2  =  u2F2<W 


where  u^  and  u-,  are  the  total  populations  of  the  two  competing  species 
and  Fj  and  F-,  are  the  exponential  growth  rates  which  depend  on  the  total 
population.  These  equations  will  be  called  the  "reaction"  equations. 

The  competitive  aspect  of  the  model  is  contained  in  the  hypothesis  that 
3Fj/9u2  and  SF^/^Uj  are  both  assumed  to  be  negative;  thus  the  growth 
rate  of  each  species  decreases  with  an  increase  in  the  competing  population. 
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B.  Example 

The  main  emphasis  Is  on  the  situation 
where  the  zero  sets  of  F^  and  F2  are 
as  pictured  in  Figure  1.  Then  the  F's 
are  negative  for  large  values  of  Uj  and 
u2  (due  to  limited  resources)  and  near 
ul  =  u2  =  0  ( due  to  the  too  low  population 

level). 

Figure  1. 

These  reaction  equations  have  seven  critical  points.  Three  of  these, 
(0,0);  (  0 ,  u2)  and  (Uj ,  0)  ,  are  attractors  and  the  point  interior  to  the 

quadrant  and  closest  to  the  origin  is  a  repeller.  So  long  as  the  points  are 
assumed  to  be  non-degenerate,  this  forces  the  remaining  three  points  to  be 
saddles . 

C.  The  Diffusion  Terms  and  Travelling  Wave  Solutions 

In  the  presence  of  diffusion,  the  model  becomes 

2  2 

auj/at  =  9  Uj/9x  +UjFj(Uj,u2) 

(2  ) 

9u2  /9t  =  62  92  u2  /9x2  +  u2  F2  ( Uj ,  u2 ) 

where  6j  and  62  are  positive  constants  measuring  the  diffusion  rates  of 
the  two  species  (these  are  not  assumed  to  be  equal). 

The  aim  is  to  show  the  existence  of  travelling  wave  solutions  for  the 
diffusion-reaction  equations.  These  are  solutions  of  the  form  u(x-0t)  . 


In  particular  they  depend  only  on  the  one  variable  £  =  x  -  0t ,  so  satisfy 
ordinary  differential  equations  of  the  form : 

-0  Uj  =  SjUj"  +u1F1(u1,u2) 

(3) 

-0U2  =  62Uz'  +  U2F2  (VU2  >  * 

These  are  most  onveniently  written  as  a  first  order  system : 


(4) 

V1  =  _6l-1  t0vl  +  uiFi  (ui»  u2^ 

v2  =  -62~:£0v2  +u2F2(ul'u2)i  * 

Travelling  wave  solutions  are  so  called  because  they  retain  their  shape 
as  time  progresses ,  but  translate  along  the  x-axis  at  velocity  8  (  since  the 

"center"  £  =  0  is  at  x  =  8t )  . 

The  waves  of  interest  are  those  which,  for  fixed  t ,  have  limiting 
values  as  x  tends  to  plus  or  minus  infinity.  This  means  that  as  solutions 
of  the  ordinary  differential  equations,  (4),  they  are  orbits  which  run  from 
one  critical  point  to  another.  The  critical  points  of  most  interest  in  this  regard 
are  the  (three)  attractors  of  the  reaction  system. 

D.  Strength  of  Attractors 

The  wave  can  be  interpreted  as  measuring  the  relative  "strength"  of 


the  attractors  as  follows :  The  system  is  in  one  stable  state  at  -®  and 
another  at  +®  ;  then  the  shape  of  the  wave  reflects  the  way  the  transition 
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from  one  state  to  the  other  will  most  naturally  take  place.  The  wave  velocity 
now  determines  which  state  will  "take  over" :  that  is,  if  ©  >  0  then  the 
state  at  minus  infinity  is  "dominant"  in  the  sense  that  the  solution  converges 
pointwise  in  x  to  this  state  as  t  tends  to  infinity.  If  0  <  0  ,  the 
other  state  is  strongest.  Of  course  the  relative  strength  defined  in  this  way 
generally  depends  on  the  way  the  diffusion  takes  place. 

E.  Relation  to  Clines 

The  existence  of  travelling  wave  solutions  depends  on  the  fact  that 
the  reaction  terms  in  the  equation  are  independent  of  x  .  However ,  the 
idea  of  the  relative  strength  of  attractors  appears  naturally  even  in  the 
x -dependent  case. 

Suppose  that  the  F' s  depend  on  x  ,  but  are  asymptotically 
space- independent  as  x  tends  to  infinity  in  either  time  direction.  Also 
suppose  the  general  appearance  of  the  zero  curves  of  F  ( for  each  fixed 
value  of  x  )  are  as  in  Figure  1  for  all  x  . 

Then  it  is  possible  that  one  attractor  dominates  at  minus  infinity  and 
the  other  at  plus  infinity,  the  dominance  being  measured  as  above  by 
the  travelling  wave  solutions  of  the  asymptotically  space  independent  systems. 
In  this  case  the  travelling  wave  is  replaced  by  a  static  solution  which  limits 
to  the  relevant  dominant  attractor  in  either  x -direction.  Such  solutions  are 
called  dines  (in  population  genetics).  This  topic  will  be  pursued  further 
someplace  else. 
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F.  Method  of  Proof 

The  relevance  of  the  remarks  about  dines  in  this  report  on  travelling 
waves  comes  from  the  method  of  proof.  Namely,  the  existence  of  either 
type  of  solution  can  be  traced  to  the  non-triviality  of  a  topological  (homotopy) 
invariant  called  the  connection  index.  The  construction  of  this  invariant  is 
described  in  detail  in  §2  (the  relevant  proofs  are  in  [1]  ,  [2]  and  [3]  ). 

The  main  job  of  the  paper  is  to  compute  this  invariant  for  the  example  in 
Fig.  1. 

The  computation  is  carried  out  much  as  analogous  computations  in 
degree  theory.  The  equations  are  deformed  or  "continued"  (  preserving  the 
invariant)  to  a  system  for  which  the  invariant  is  easy  to  compute.  This  simpler 
system,  described  in  §3,  is  itself  the  product  of  a  "hyperbolic  point"  and  a 
lower  dimensional  (planar)  system.  A  product  theorem  for  the  invariant  thus 
reduces  the  computation  to  that  for  a  planar  system. 

The  planar  system  is  described  in  §  2  ( with  some  detail  deferred 
to  the  appendix) ;  it  is  also  used  there  to  illustrate  the  abstract  definition  of 
the  connection  index. 

In  §4  "isolating  neighborhoods"  (analogous  to  open  sets  with  no  fixed 
points  on  the  boundary)  are  constructed ;  these  are  needed  to  carry  out  the 
continuation  described  in  §5.  This  construction  illustrates  some  novel  ways 
of  locating  Isolating  neighborhoods  in  higher  dimensional  systems  of  this  type. 
The  lemmas  proved  are  somewhat  more  general  than  necessary ,  but  allow  for 
applications  to  examples  other  than  that  of  Figure  1. 

G.  An  Alternate  Proof  and  Stability 

In  a  forthcoming  paper  by  the  second  author  ([4])  a  proof  of  existence 
of  the  travelling  wave  will  be  obtained  from  a  topological  degree  argument. 
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It  seems  clear  that  the  success  of  one  approach  implies  that  of  the  other 
in  some  general  setting;  however,  a  general  theorem  has  not  been  formulated. 

A  similar  situation  occurs  in  the  study  of  shock  waves  (compare  [5]  and  [6]). 

Also  in  the  above  mentioned  paper,  [4],  the  stability  of  the  waves 
as  solutions  of  the  partial  differential  equation  will  be  treated. 

H.  The  Main  Theorem. 

In  the  example  of  Figure  1 ,  there  are  three  attractors  :  namely 
( u2 , 0)  ,  (0,0)  and  ( 0  ,  u^)  .  It  is  already  well  known  ( and  the  same  invariant 
has  been  used  in  [2]  to  see )  that  there  are  travelling  waves  running  from  the 
first  to  the  second  and  from  the  second  to  the  third  (these  are,  respectively, 
solutions  of  the  two  dimensional  systems  obtained  on  setting  Uj  or  u2 
equal  to  zero) . 

The  main  theorem  concerns  waves  from  (u2,0)  to  (0,1^)  : 

Theorem :  With  reference  to  A.,  B.  and  the  above  remarks ,  let  0' 

and  0"  be  the  wave  velocities  of  the  waves  from  (u2 ,0)  to  (0,0)  and 
from  (0,0)  to  (O.up  respectively.  Then  if  0'  >  0"  ,  there  is  also  a 
wave  from  (u2,0)  to  ( 0 ,  Uj)  . 

The  proof  of  this  theorem  is  concluded  in  §  5  . 
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§  2  .  The  Standard  Problem; 

As  described  in  the  introduction,  the  connection  problem  treated  here 
can  be  continued  to  a  "standard"  problem  which  is  described  below.  This 
model  problem  will  also  be  used  to  illustrate  the  abstract  definitions  required 
to  arrive  at  the  connection  index  used  in  the  existence  proof. 

The  equations  are 


(5) 


w '  =  z 

2 

z1  =  -8z  -  w(l-w  )  . 


The  phase  portrait  in  the  (w,  z) -plane  with 
Figure  2 . 

There  is  a  solution  in  the  upper  half 
plane  connecting  the  left  hand ,  hyperbolic 
rest  point,  p'  =  (-1, 0)  ,  to  the  right  hand 
one ,  p' 1  =  ( 0 , 1)  .  A  homotopy  invariant 
related  to  the  existence  of  this  connection 
is  constructed  as  follows. 


9=0  is  sketched  in 


*  * 


Figure  2 

! 

* 

A.  Definition  1. 

Suppose  a  differential  equation  (or  local  flow)  on  a  space  X  is 
given.  If  N  C  X  is  compact,  let  I(N)  denote  the  set  of  points  on 
solutions  which  stay  in  N  for  all  time.  (Then  I(N)  is  a  compact 
invariant  set). 
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If  S  =  I(  N  )  is  interior  to  N  ,  it  will  be  called  an  isolated  invariant 
set  and  N  will  be  called  an  isolating  neighborhood  for  S  . 

Definition  2 . 

Now  suppose  there  is  given  a  family  of  differential  equations  on  Rn 
parametrized  (continuously)  by  a  parameter  9  in  a  closed  interval  [6^0^] 
Let  X=  Rn  x  £©0 ,  ©j]  .  Then,  in  view  of  the  continuity,  the 
equations  determine  a  (local)  flow  on  X  . 

Let  S' ,  S"  and  S  be  isolated  invariant  sets  for  this  local  flow 
and  let  S'  ( 0)  etc.  be  the  set  of  points  in  S’  with  parameter  value  9  . 

Then  the  triple  (  S' ,  S"  ,  S)  will  be  called  a  (  codimension  1 ) 
connection  triple  if 

a)  S'  U  S"  C  S 

b)  S'  n  S"  =  9 

c)  For  0  =  0Q  or  9^  S(0)=  S' ( 9)  U  S' ' ( 9)  . 


B.  Example . 


With  regard  to  Figure  2  ,  let  D'  and  D’ 1  be  disjoint  compact 
neighborhoods  of  p'  and  p"  respectively,  neither  of  which  contains  (0,0) 
Let  D  be  a  compact  neighborhood  of  the  connecting  solution  disjoint  from 
the  non-positive  z-axis.  Then  D'  and  D"  are  isolating  neighborhoods 
for  all  values  of  9  with  corresponding  isolated  invariant  sets  p'  and  p' '  . 

Also,  D  is  an  isolating  neighborhood  for  all  0  :  if  0  /  0  , 

1(D)  =  {  p'  ,  p"}  and  if  0=0,  1(D)  is  the  closure  of  the  connecting 

solution. 

It  follows  that  for  any  compact  interval  [0Q,  0^]  ,  N'  s  D'  x  [0q,0^]  , 
N"  =  D"  x  [6q*®i]  and  N  =  Dx  [0q,0^]  are  all  Isolating  neighborhoods. 
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Also  S'  =  I(N')  and  S"  =  I(N")  are  disjoint  and  contained  in  S=  I(N)  , 
so  if  condition  c)  of  Definition  2  is  satisfied  then  (  S' ,  S"  ,  S  )  is  a 
connection  triple.  This  is  the  case  if  neither  of  0Q  nor  0^  is  0  ; 
however,  if  0Q  (say)  is  0  then  S( 0Q)  5^  S' ( 0Q)  U  S"(0q)  and 
condition  c)  is  violated. 

Observe  that  if  0Q  and  0^  have  the  same  sign  then  S  =  S’  U  S” 
and  there  is  no  connection  (for  any  0  in  [0q,9^])  from  S’  to  S"  . 
However,  if  0Q  and  0j  have  different  signs  then  S  /  S’  U  S"  and  the 
difference  S\S'U  S"  is  the  connecting  solution.  The  aim  is  to  define  a 
homotopy  invariant  for  a  connection  triple  which  will  distinguish  these  two 
cases.  Before  defining  this  invariant,  '’continuation”  (or  deformation)  of 
isolated  invariant  sets  and  triples  will  be  defined. 

C.  Definition  1 . 

Suppose  now  that  there  is  given  a  family  of  local  flows  on  a  space  X 
parametrized  (continuously)  by  a  parameter  X  t  [0,1]  .  Let  Y  =  Xx  [0,1] 
and  consider  the  family  as  a  local  flow  on  Y  .  Given  a  compact  set  N  C  y 
and  X  «  [0,1]  ,  let  N  (X)  be  the  set  of  points  in  N  with  parameter  value  X 
( Observe  that  if  N  is  an  isolating  neighborhood ,  then  for  X  e  [0,1]  , 

N  (X)  is  also). 

Let  Sq  and  be  Isolated  invariant  sets  for  the  flows  on  X 

corresponding  to  X  =  0  and  X  =  1  respectively. 

Then  SQ  and  Sj  are  related  by  continuation  if  there  is  an  isolating 
neighborhood ,  N  ,  for  the  flow  on  Y  such  that  SQ  =  I  (  N  ( 0) )  and 
Si  =  I(Nj)  .  (One  could  say  that  as  X  varies  from  0  to  1  ,  SQ  is 
"deformed"  to  Si  .  However,  it  is  really  the  equations  that  are  deformed— 
the  sets  may  change  radically). 
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Definition  2 

In  a  like  way  one  can  define  continuation  of  connection  triples : 
assume  that  in  the  above,  X  =  Rn  x  [0Q ,  0^]  and  that  for  each  X  ,  the 
flow  on  X  is  actually  a  family  of  flows  parametrized  by  0  .  Suppose 
also  that  S  ,  S'  ,  S"  are  isolated  invariant  sets  of  the  flow  on  Y  and 
that  for  each  X,  (S' (A),  S"(\) ,  S(  A.) )  is  a  connection  triple.  Then 
the  triples  at  X  =  0  and  X  =  1  are  related  by  continuation. 

The  " continuation"  relation  is  obviously  an  equivalence  relation; 
the  aim  now  is  to  define  an  invariant  or  "index”  which  is  the  same  for  all 
isolated  invariant  sets  (or  connection  triples)  that  are  in  the  same  equivalence 
class. 

D.  Definition. 

Let  N  be  an  isolating  neighborhood  for  a  flow  on  a  space  X  and 
let  Nj  and  N2  C  be  compact  subsets  of  N  .  Suppose  Nj\N2  is 
a  neighborhood  of  I  ( N)  and  that  Nj  and  N2  are  positively  invariant 
relative  to  N  .  (I.e.  whenever  the  initial  point  of  an  orbit  segment  lies  in 

Nj  ,  respectively  N-,  ,  and  the  segment  lies  in  N  ,  then  the  segment  lies 
in  Nj  ,  respectively  N2  ).  Also  assume  that  points  of  which  are 

carried  out  of  N  hit  N2  before  leaving  N  . 

Then  (Nj,N2)  is  called  an  index  pair  for  S=I(N). 

Theorem  1.  For  any  isolated  invariant  set  S  ,  index  pairs  exist.  Moreover 
the  homotopy  type,  [Nj,N2]  ,  of  the  pointed  space  N|/ N2  depends  only  on 
S  =  I  ( N )  and  is  the  same  for  all  isolated  invariant  sets  which  are  related  by 
continuation  to  S  .  It  will  be  called  the  (homotopy)  index  of  S  and  denoted  h(S) 


fc  i -ai- 
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Remark :  If  SQ  and  are  related  by  continuation  ( as  In  Definition  1 

of  C.  )  and  N  is  the  neighborhood  used  to  define  the  continuation  then 
h(SQ)=  h  ( I  ( N) )  . 

The  definition  of  the  connection  index  is  played  back  to  that  above 
using  the  following  lemma. 

Lemma  :  Let  ( S' ,  S“  »  S )  be  a  connection  triple  for  a  family  of  differential 

equations  on  Rn  parametrized  by  0  in  the  interval  [0Q,0j]  •  Assume 
the  equations  are  defined  for  0  c  (0^  -  e  ,  0j+e)  for  some  e  >  0 
(this  is  no  real  restriction— they  can  be  extended  to  such  an  interval).  Let 
U‘  and  U"  be  open  neighborhoods  in  Rn  x  (0q-e,  0^+e)  of 
S'(0q)  U  S'(0^)  and  S"(0q)  U  S'^Oj)  (respectively);  choose  these  to 
have  disjoint  closures.  Let  <p  be  a  continuous  real  valued  function  on  Rn 
which  is  positive  on  U  1  and  negative  on  U  ”  .  Append  to  the  given  family 
of  equations  the  equation  0'  =  p?>(x)  (0Q  +  0j)/2)  where  \i  is  a  (small) 
positive  parameter. 

Then  there  is  a  nQ  >  0  such  that  if  \i  e  (0,^Q)  then  N  is  an 
isolating  neighborhood  for  the  appended  equation.  Let  h^  be  the  index  of 
S(N)  ,  p  «  (0,n0)  .  Then  h^  is  independent  of  p  ,  and  in  fact 
depends  only  on  the  triple  ( S' ,  S"  ,  S  )  . 

Proof.  This  lemma  is  proved  in  [3]  in  a  rather  more  general  form  ( see  also 

[2])- 

Definition :  With  the  above  notation ,  the  "  index"  of  the  triple  ( S' ,  S"  ,  S ) 
is  h(S',S",S)=h  . 
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Theorem :  The  index  of  a  connection  triple  is  constant  on  equivalence 

classes  under  the  continuation  relation. 

Proof.  See  [3]  . 

E.  Example.  ( See  Appendix  A ) 

The  indices  of  the  hyperbolic  rest  points  in  (5)  (Example  B.)  are 
both  pointed  circles,  denoted  here  by  E1  . 

The  index  of  the  disjoint  union  of  isolated  invariant  sets  is  the  "sum" 
of  their  indices;  namely,  the  pointed  space  obtained  on  identifying  their 
distinguished  points.  Thus  the  index  of  the  union  of  the  two  hyperbolic 
points  is  E1  V  E1  where  V  means  sum. 

By  the  continuation  theorem,  the  index  of  the  closure  of  the  connecting 
solution  is  also  E1  v  E1  .  (This  set  is  the  continuation  of  the  set  {  p' ,  p"} 
with  e  e  [0Q,O]  the  continuation  parameter). 

The  index  for  the  connection  triple  in  B.  is  shown  in  Appendix  A 
to  be  E2  v  E*  if  0Q  and  0^  have  the  same  sign  (i.e.  when  0=0 
and  the  connecting  solution  are  absent )  and  TT  ( the  homotopy  type  of  the 
one  point  space)  when  0Q0^  <  0  . 

Remark ;  In  addition  to  the  sum  of  indices,  a  product  is  also  defined  by 
[Nj/N^]  A  [Nj"  ,  N^']  [Nj'xNj/iNjxNj'  l)  N^'xNj}]  .  Given  a  product 
of  flows ,  the  product  of  isolated  invariant  sets  is  isolated  and  its  Index  is  the 
product  of  the  Indices  of  the  sets.  This  product  is  used  in  the  statement  of 
the  "existence"  theorem  following. 


F.  The  theorem  needed  for  the  existence  proof  is : 


Theorem :  Suppose  (S' ,  S"  ,  S)  is  a  connection  triple  with  index  h  , 

and  let  h'  and  h"  denote  the  indices  of  S'  and  S’ 1  respectively. 
Then  if  h/ (Z1  A  h’)  V  h"  ,  S/S'US"  . 

Example :  (Seethe  appendix).  In  the  example  (E.)  if  0^9^  >  0 

then  h  =  Z2  v  Z1 ,  h'  =  h"  =  Z1  ,  and  Z1  A  Z1  =  Z2  .  Therefore 
h  =  (  2*  A  h')  v  h"  so  the  theorem  gives  no  information.  If  9O01  <  0 
then  h  =  0  /  S2  v  2^  so  S/S'  U  S"  .  The  difference  between  S  and 
S'  U  S"  is,  of  course,  the  connecting  solution  in  this  case. 

G.  In  §4  and  §  5  it  will  be  seen  that  the  problem  of  §  1  can  be  resolved 

by  continuing  an  appropriate  connection  triple  to  a  product  of  that  in  B.  with 

a  hyperbolic  rest  point.  The  indices  of  the  triple  in  the  product  will  then  be 

h =  Z1 A  0  =  0  ,  h'  =  Zl  A  Z1  =  Z2  and  h"  =  Z2  .  Then 
“  1  2  2  3  2 

0/  (Z  A  Z  )VZ  ( =  Z  VZ)  will  imply  the  existence  of  the  connecting 
solution. 

The  continuation  relies  only  on  identifying  isolating  neighborhoods 
which  is  fairly  easy.  Before  doing  this,  however,  the  "target"  equation 
(i.e.  the  product)  will  be  described  in  §3. 


ik 


i  3 

V1  =  "0V1  “  ^4ul  “  6u1u2  ~  4ul  ) 

i  3 

v2  =  “0V2  '  (4u2  "  6U1U2  -  4U2>  * 

The  zero  sets  of  Fj  and  F2  are  pictured  in  Figure  3 


Figure  3, 


Wj  =  Uj  +  u2  and  w2  =  U1  ”  u2  ’  Then 

2  2  2V  2  .  2 
w2  ,  2  { Uj  +  u2  )  =  Wj  +  w2  , 

W1  l2(wi  +  w2)  -  (v/2  -  w2)]  =  Wj  [w2  +  3w2  ]  i 
w2  [2(w2  +  w2)  +  (w2  -  w2))  =  w2[3w2  +  w2]  . 


Let 

4u1u2  =  - 

4(u13+  u2)  = 

.  .  3  3, 

4  ( U|  -  u2)  = 


The  transformed  equations  are  therefore 
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W1  =  Z1 
Z2 

-02j  -  4wj  +  3(w2  -  w 2 )  +  Wj(w2  +  3w2  ) 

•  2  2 
z2  =  -0z2  -  4w2  +  w2(3wj  +  w2)  . 


(7) 


W-,  = 


z,  = 


B.  An  Invariant  Manifold  Containing  the  Standard  Problem. 

When  Wj  =  1  and  z^  =  0  ,  wi  =  z,1  =  0  .  Therefore  the  set 
{wj  =  1 ;  Zj  s  0}  Is  an  invariant  manifold.  On  this  manifold,  the  equations 
for  w2  and  z2  are : 


(8) 


w2  3  z2 


*  2 
z7'  =  -0Z7  -  w7(l  -  w7)  ; 


namely  the  equations  of  §  2  . 

The  displacement  equations  describing  the  behavior  of  Wj  and  z^ 
near  the  invariant  manifold  are  given  by : 


(9) 


i  =  ^ 


=  -0;  +  (5  +  3w2)S  +  0(42) 


Wj  -  1 1  2=  e/K 


(wherein  Wj  =  1  +  4  and  Zj  =  4  )• 

Suppose  now  that  N  is  any  compact  set  in  ( w2 ,  z2 , 0)  -space  . 

I 

Let  K  be  a  bound  for  6w2  z2  on  N  .  Let  D£  =  {(w1#z^) 
and  |zj|  s  e  .  The  idea  is  to  show  that  if  e  is  small  enough,  then 
I  (D£  x  N)  =  {(1,0)}  x  I ( N)  . 

To  see  this  it  is  sufficient  to  see  that  for  small  e  ,  no  solution  in 
I(D  x  N)  can  pass  through  a  point  in  00  x  N  . 

E  £ 


Suppose,  then,  that  |4|  =  e/K  .  If  4*  =  4  is  not  zero  the 

solution  can't  be  in  I  (D  x  N )  .  If  4  =  0  »  then 

4"  =  -04  +  (  5  +  3w^ )  4  +  0  (42)  •  Thus  4"  has  the  same  sign  as  4> 

so  again  the  solution  is  not  in  I(D£  x  N)  if  e  is  small  enough. 

Now  suppose  |4|  =  e.  If  4*  =  4  /  0  then  the  solution  isn’t  in 

I(D  x  N)  .  If  4=0  then  4"  =  4  so  4"  has  the  same  sign  as  4 

and  again  the  solution  must  leave  Dg  x  N  . 

Observe  that  the  same  argument  as  that  above  works  if  the  terms  in 
2  2 

( 9 )  involving  0  ,  w2  and  4  are  multiplied  by  a  parameter  X  varying 

between  0  and  1  . 

Now  let  N  be  an  isolating  for  the  ( w2  ,  z2 , 0)  -equation. 

Then  for  small  e  ,  D£  x  N  is  an  isolating  neighborhood  for  (  6)  or 
equivalently,  (8)  with  (9)  .  Namely, 

I  ( D£  x  N  )  =  I  (  {  0}  x  N  )  C  int.  D£  X  N  .  Modifying  the  equations  with 
a  parameter  X  as  described  above,  one  continues  the  invariant  set  to 
(the  "same")  one  for  the  product  of  the  equation  (8)  with  the  hyperbolic 
point : 


(10) 

41  =54  (x  =  0) 

Thus  the  connection  triple  in  the  invariant  manifold  Wj  =  1  ,  =  0 

with  0  c  [0Q,  0j]  continues  to  the  product  of  that  in  2.B.  with  a 
hyperbolic  point.  As  described  in  2.G.,  when  0^1  c  0,  the  indices 
for  triple  are  h  =  0  ,  h'  =  22  and  h"  =  22  . 

In  the  next  section,  isolating  neighborhoods  will  be  identified  which 
allow  one  to  continue  the  general  problem  to  the  one  treated  here. 
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§  4.  The  Main  Problem 

The  equations  of  §  1  are  written  here  In  the  more  convenient  form : 

U1  =  V1 

U2  =  v2 

(ID 

V1  =  -91V1  "  ulfl (  ul*  u2> 

V2  =  “02V2  "  u2f2(u  ' 

where  9^  =  6”1 9  and  fj  =  6”1  Fj  (1=1,2). 

Observe  that  9^  and  0 2  have  the  same  sign  (or  are  zero  together); 
this  is  necessary  for  the  results  on  isolating  neighborhoods  that  follow. 

Also  note  that  and  6-,  do  not  appear  explicitly  in  the  following ; 
the  aim  is  to  continue  a  connection  triple  of  the  equations  (11)  above  to  that 
for  equations  (6)  of  §3,  and  in  so  doing  it  does  not  matter  that  the  relation 
between  the  6^  and  0^  is  lost. 

To  carry  out  the  continuation,  it  is  sufficient  to  find  appropriate  isolating 
neighborhoods,  or  in  other  words,  compact  sets  N  such  that  I  ( N)  n  9N  =  <P  . 
This  will  be  done  by  "successive  approximations." 

The  following  hypotheses  hold  throughout; 

Hypotheses.  There  is  a  positive  constant  K  such  that  if  u  lies  in  the 
quadrant  u4  s  0  (1=1,2),  then: 

1)  Uj  2  K  implies  fj(Uj,u2)  <  0  ,  and  u2  >  K  implies 
f2  ( u^ ,  u2)  <  0  . 
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2)  There  are  only  finitely  many  points  u  for  which  f^  =  f 2  =  0  . 

3)  3f1/au2  and  3f2/3uj  are  both  negative. 

A.  Definition . 

Let  S  be  the  set  of  points  on  bounded  solutions  of  ( 11)  whose 
u -coordinates  always  satisfy  Uj  a  0  (1=1,2). 

Lemma :  There  is  a  constant  L  >  0  such  that  for  any  values  of  0^  and 

02,  SC  {(u,v)|  0  s  Uj<  K  and  IvJ  <  L  1=1,2}  . 

Proof.  Choose  any  solution  in  S  and  let  ( u ,  v)  be  a  point  in  its 
closure  where  Uj  achieves  its  maximum  value.  Then  Uj  =  Vj  must  be 
zero  so  uj'  =  -Ujfj(Uj,u2)  .  But  also,  uj'  must  be  non-positive  so  fj 

must  be  non-negative.  This  implies  u^  <  K  ,  and  it  follows  that 
Uj|  S  <  K  .  A  similar  argument  proves  u2  |  S  <  K  . 

Now  let  M  be  a  bound  for  |uifi  (u^,  u2)  |  on  the  square 
0  s  Uj  s  K  (i  =  1,2)  .  Suppose  a:  1  .  Consider  a  solution  on  which 
|  vj  exceeds  M  at  some  time.  Then  at  that  time  Vj  and  v^  have  the 
same  sign  and  this  condition  will  be  maintained  either  in  forward  time 
( Vj  >  0  )  or  backward  time  (v^  <  0)  so  that  the  solution  will  grow 
exponentially  and  not  be  in  S  .  Similar  arguments  show  that  if  joj  a  l 
then,  on  S  ,  |vj|  s  M  . 

Suppose  |0j)  £  1  .  Then  for  large  enough  Vj  ,  |vj  |  <  2|vj|  . 
Therefore,  there  is  a  constant  L'  >  K  such  that  if  |vj  is  larger  than  L' 
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at  time  zero,  it  is  larger  than  K  for  time  1  .  It  follows  that  at  zero 
differs  from  Uj  at  1  by  more  than  K  and  the  solution  cannot  be  in  S  . 
Thus  on  S  ,  |vj|  ,  and  in  a  similar  way  ( v2  |  ,  is  strictly  bounded  by 
L  =  max  (  M  +  l,  L'  )  independently  of  8j  and  02  • 

B.  Definition . 

Let  C  denote  the  set  of  critical  points  in  S  .  Let  Uj+  denote 
the  points  on  solutions  in  S  with  u-,  =  v2  =  0  and  Vj  2  0  and  U2 
those  with  Uj  =  =  0  and  v2  s  0  . 

Lemma  :  Let  s  { (u ,  v) j  0  s  u^.^sK  and  0  s  ,  - v 2  ^  L }  . 

Then  I  ( Nj)  consists  of  critical  points  and  orbits  connecting  them. 
Furthermore,  if  p'  «  9Nj  then  either 

a)  p'  is  on  an  unbounded  solution,  or 

b)  The  solution  through  p'  leaves  Nj  immediately  in  one  or 
the  other  time  direction ,  or 

c)  The  projection  in  u- space  of  the  solution  through  p'  is 
contained  in  one  of  the  coordinate  axis,  or 

d)  p'  «  C  . 

In  particular ,  I  ( Nj)  n  9Nj  =  U  uj"  U  C  . 

Proof.  Since  Vj  2  0  and  v2  s  0  in  Nj  ,  the  function  Uj  -  u2  is 

strictly  increasing  on  non-constant  solutions  in  .  (If  Uj  -  u2  were 
constant  on  some  time  interval,  both  v^  (  a  0  )  and  v2  (  s  o  )  would  be 

zero  on  the  interval  and  the  solution  would  be  a  critical  point).  By  general 


theorems,  it  follows  that  I(N^)  consists  of  critical  points  and  solutions 
connecting  them  (cf.  [1]). 

Suppose  p'  e  3Nj  .  Then  either  some  u^  is  K  or  0  ,  or  one 
of  Vj  or  -Vj,  is  L  or  0  . 

1.  If  Uj  or  u-,  equals  K  ,  or  v^  or  -v^  equals  L  ,  the  solution 

through  p'  is  unbounded  {by  A.  ). 

2 .  Suppose  Uj  =  0  .  If  uj  =  is  not  zero  the  solution  leaves 

immediately  in  one  or  the  other  time  direction. 

If  Vj  =  0  the  projection  of  the  solution  in  u -space  is  in  the  u2 
axis  (since  the  set  u^  =  v^  =  0  is  obviously  an  invariant  set). 

The  case  u^  =  0  is  treated  similarly. 

3.  Suppose  Vj  =  0  and  Uj  /  0  .  If  Vj  =  -u^fj  Is  not  zero  the  solution 

leaves  immediately  in  one  or  the  other  time  direction.  If  v^  =  0 

then  v^  =  -Uj  ( 9fj/9u2)  v2  .  If  this  is  not  zero  then  v^"  <  0 

(Uj  >  0  ,  df^/du^  <  0  and  v^  <  0)  so  v^  has  achieved  a  strict 
maximum  and  the  solution  leaves  Nj  immediately  in  both  time  directions. 

4.  Continuing  from  3. ,  suppose  v^'  =  0  ;  then  v2  =  0  .  If  uz  =  0 

the  projection  is  in  the  u^  axis.  Assume  u.,  /  0  .  If  v^  /  0  the 

solution  leaves  Ny  immediately  (as  in  3.  ).  Otherwise 

Vj  =  v,  =  v^  =  v^  =  0  and  the  point  is  in  C  • 

This  concludes  the  proof  of  Lemma  B. 


C.  Increasing  Near  Attractors  of  the  Reaction  Equations. 

Let  p=  (Uj,u->,0,0)  be  a  critical  point  of  (8)  and  let  Mp  be 
the  Jacobean  matrix  9(Uj  fj,  u2  f2)/3(  Uj ,  u2)  at  p.  Let  Cj  be  the  set 
of  critical  points  p  such  that  the  diagonal  entries  of  Mp  are  negative 
and  det  Mp  >  0  .  (These  points  are  easily  seen  to  be  attractors  for  the 
reaction  equations  (1)  of  the  introduction ,  but  that  fact  will  not  be  used 
explicitly) . 

Given  positive  numbers  Aj  and  A^  and  p  «  [0,1]  ,  define 
N  (p,  P)  s  {(U,  v)|  |  Uj  -  Ujl  Aj  ;  |u2  -  u2|  ^  A2  and  |  v±  |  <  L  ,  i  =  1 ,  Z  > 

Lemma :  Positive  numbers  A^  and  A^  can  be  chosen  so  that  for  all 

p  c  (0,1]  ,  N(p,p)  is  an  isolating  neighborhood  for  p  . 


Having  chosen  numbers  (A^.A^)  as  above  for  each  point  of  Cj. 
so  that  the  sets  N(l,p)  (p  e  Cj)  are  disjoint,  let  N2(p)  be  the  union 
of  Nj  with  all  the  N(p,p)  for  p  «  Cj  .  Let  N->  =  N2(l)  . 

Then  for  p  e  (0,1]  ,  I(N2)  =  I(Nj)  •  In  particular,  the  points  of 

Cj  are  not  in  I(N2)  D  9N2  . 

Proof. 

1.  Since  the  off-diagonal  terms  of  Mp  are  known  to  be  non- positive ,  the 
hypotheses  imply  that  the  gradient  vectors  of  Uj  fj  and  u^  at 

(Uj,  u^)  both  point  into  the  closed  third  quadrant  and  that  the  orientation 
determined  by  the  ordering  (V(U|fj),  V u^  b>)  agrees  with  the  standard 
orientation  of  the  (u^,  u2)  plane. 

2.  In  view  of  1.  ,  there  is  a  line  with  negative  slope  through  (u^,  u2)  such 
that  at  points  on  the  line  near  (u^,u2)  ,  both  of  ( u^ -  u^) u^  f^  (u^,  u2) 
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and  (u2  -  u2)u2  f2(Uj,  u2)  are  non-positive  and  if  either  one  is  Eero 
then  u^  =  Uj  and  u2  =  u2  .  (See  Figure  4). 


■J.  U  =  o 


Figure  4 

Let  this  line  be  given  by  A^  ( u^  -  u^)  +  A2(u2  -  u2)  =  0  where  Aj 
and  A2  are  positive  and  chosen  small  enough  that  the  above  inequalities 
are  satisfied  on  the  segment  of  the  line  determined  by  |  u^  -  uj  £  A^ 
(1=1,2)  (Fig.  4). 

3.  Now  choose  ^  t  (0,1]  .  The  u -coordinates  of  points  in  N(|i,p) 

then  range  over  the  rectangle  |u^-u^|  s  u  A^  ;  |  u2  -  u2  |  ^  fiA2  . 

Observe  that  u^  fj  is  positive  on  the  left  hand  side  of  the  rectangle  and 
negative  on  the  right,  while  u?  f2  is  negative  on  the  top  and  positive 
on  the  bottom. 

4.  To  see  that  N  (n,  p)  is  an  isolating  neighborhood,  first  recall  that  no 
orbits  on  which  Vj  or  v2  equals  L  are  bounded.  Therefore  any 
point  in  I  (  N  (  n ,  p) )  n  9N  (p. ,  p)  must  have  u  -coordinates  in  the 
boundary  of  the  rectangle-  Also  the  orbit  through  the  point  must  be  tangent 
to  the  boundary  (cusp  points  counting  as  tangencies). 
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Suppose ,  for  example ,  that  the  point  is  in  the  left  hand  boundary 
so  that  u^  —  TT^  =  -p  and  =  0  .  Then  v^'  =  -u^f^  and  is  negative. 

It  follows  that  -  u^  has  a  strict  maximum  at  the  point  and  from  there 
that  the  solution  leaves  N  ( pt)  in  both  time  directions  so  in  particular 
the  point  can't  be  in  I  (  N  (p ,  p) )  fi  3N  (p  ,  p)  . 

The  other  three  cases  are  treated  in  a  similar  way  and  it  is  found  that 
if  |u^  -  u,|  =  pA^  for  i  =  1  or  2,  then  the  orbit  through  u  leaves 
N  (p ,  p)  in  both  time  directions. 

5.  Now,  since  N(p,p)  is  an  isolating  neighborhood  for  all  p  €  (0,1]  , 
the  corresponding  isolated  invariant  set  must  be  in  the  intersection  of  the 
N  ( p ,  p)  .  The  points  in  this  intersection  have  the  form  ( Uj ,  u2  ,  Vj ,  v2) 
so  the  only  orbit  in  the  intersection  is  the  critical  point. 

6.  To  prove  the  last  paragraph  of  the  lemma,  suppose  N2(p)  is  defined  as  it  is 
there  and  choose  p  so  that  the  N(p,p)  (p  e  Cj)  are  disjoint.  It 

will  be  shown  that  I(N2(p))  H  9N2(p)  C  I(N^)  for  all  p  . 

Assuming  this  for  now,  one  sees  that  I(Nj)  =  I(N2(p))  as  follows. 
First  let  N2(0)  =  O  (N2(p)i  p  «  (0,1]}  .  Then  points  in  N2(0)\Nj 
have  the  form  (u^,  u-, ,  v^,  v2)  where  (u^,u2,0,0)=  p  <•  Cj.  and  v 
violates  the  condition  0  ==  v^  ,  -v2  £  L  .  These  points  are  obviously 

not  in  I  (N2(  0) )  (  since  they  leave  N2  (  0)  immediately  in  both  time 

direction)  and  so  I(N2(0))  =  I(Nj)  . 

Now  if  I  (N2(p) )  /  I  (Nj)  for  some  p  ,  then  for  some  p  , 

I(N2(p))  n  9N2(p)  contains  a  point  not  in  I(Nj)  .  But  this  violates 
I(N2(P))  n  9N2(p)  c  KNj)  . 

7.  To  show  that  for  all  p«  (0,1],  I(N2(p|)  fl  3N2(p)  C  I(Nj)  ,  first 

observe  that  9N2(p)  consist  of  three  types  of  points:  Namely, 
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a)  points  not  in  Nj  but  in  the  boundary  of  some  N(p,p)  with  p  i  Cj  ; 

b)  points  in  3N^  but  not  in  any  N  ( p ,  p)  ;  and 

c)  points  in  3Nj  n  3N  ( p ,  p)  for  some  p  «  Cj  . 

8.  Suppose  p'  c  3N(p,p)\Nj  for  some  p  c  Cj  .  Since  the  solution  leaves 
N  (p,  p)  immediately  in  one  or  the  other  time  direction,  it  leaves  N->(p) 

(it  can't  enter  another  N  (p ,  p)  since  they  are  disjoint). 

9.  Suppose  p'  e  3Nj  but  not  in  any  of  the  N(p,p),  p  t  Cj  .  Using 
Lemma  B,  it  is  either  on  an  unbounded  solution  (  so  leaves  N2(p)  )  or 
leaves  Nj  immediately  in  some  direction  (again  leaving  N2(p)  also 
since  it  is  in  dNz  (  p)  ;  cf.  4.)  or  it  is  a  critical  point  (  so  in  I  (N^) ) 
or  the  u -projection  of  the  solution  through  p'  lies  in  one  or  the  other 
of  the  coordinate  axis . 

Suppose  for  example  that  it  lies  in  the  u^  -axis .  If  v^  is  never 
negative  on  the  solution,  it  is  in  I(N^)  . 

On  the  other  hand  if  v^  is  sometimes  negative  then  (since  it  is  also 
sometimes  positive)  the  solution  must  be  an  oval,  or  spiral  down  to  a  rest 
point  in  one  or  the  other  time  direction  (cf.  Figure  5).  Of  course  all  the 
spiraling  must  take  place  near 
some  critical  point,  p  ,  in  the 
Uj  axis  since  v^  is  not 
negative  in  .  But  solutions 
in  the  U|  axis  near  p  «  Cj 

do  not  spiral,  but  rather  look  F,€i  * 

locally  like  hyperbolas  (cf.  Fig.  5). 
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10.  To  finish  the  proof  It  must  be  shown  that  the  points  In 

I(N2(p))  n  ONj  fl  N(p,p)  {for  some  p  c  Cj)  are  also  in  I { N^)  . 

In  9-  it  was  seen  that  points  in  I(N2)  with  or  u2  equal 

to  zero  are  also  in  I(N^)  .  Therefore  assume  Uj,u2  >  0  • 

Let  p=  (u^,u2,0,0)  be  the  critical  point.  If  p1  e  9Nj  n  9N(p,p) 

then  one  of  the  v^  is  zero  and  one  of  the  |uj  -  uj  is  pA^  .  If, 

Vj  =  0  and  { u^  —  u^ j  =  pA^  then  the  point  leaves  and  N(p,p) 

in  both  time  directions  (by  4.)  so  the  point  isn't  in  N2(p)  .  Similarly 
if  v2  =  0  and  ju2  -  u2|  =  pA^  ,  the  orbit  through  p'  is  not  in  N2(p)  . 

11.  Suppose  then  that  Uj  -  u^  =  -pAj  and  /  0  .  Then  v2  =  0  . 

There  are  two  cases  (indicated  by  the  curved  arcs  in  Figure  4).  Either 

u2  f2  ~  0  ^the  point  lies  on  or  bel°w  the  curve  u2  f 2  =  0  )  or  u2  f2  <  0  . 

In  the  first  case  v2  becomes  positive  as  time  decreases,  so  the 
point  leaves  as  well  as  N  (p,p)  in  negative  time.  (If  u2  f2  =  0  , 

v2'  =  -u2  9f2/9u^v^  <  0  ). 

Suppose  then  that  u2  f2  <  0  .  Then  v2  becomes  positive  as  time 
increases.  Now  following  the  solution  forward  in  time,  it  enters  N  (p,  p) 
and  u2  increases  initially.  Since  the  curve  u2  £,  has  non- positive  slope 
in  N  (p,  p)  ,  the  solution  stays  above  this  curve  at  least  until  v2  becomes 
zero.  Suppose  v2  does  become  zero  at  some  point  on  the  solution  curve. 
Choose  the  first  such  point.  At  this  point,  v2  =  -u2  f-,  is  positive.  But 
this  contradicts  the  fact  that  v2  is  positive  up  to  this  point. 

Therefore  v2  cannot  become  zero  so  long  as  the  curve  is  in  N(p,p)  . 
This  implies  u2  is  increasing  on  the  curve— so  it  cannot  approach  the 
critical  point  p  ,  and  from  there  that  the  curve  must  sometime  leave 
N(p,  p)  .  But  then  it  leaves  with  v2  positive  so  it  also  leaves 
Nl  U  N(p,p)  . 


26 


12.  The  other  cases  (uj  -  =  n  Aj ,  Vj  >  0  and  It^  -  u2|  =  , 

v2  <  0  )  are  similarly  treated ,  and  It  is  shown  that  points  In 
9N(|i,p)  fi  9N^  with  Uj,u2  >0  are  not  contained  in  N2(p)  . 

Together  with  8.  and  9.  this  completes  the  proof  that 
I(N2(h))  n  9N2(ji)  C  I  ( Nj)  ,  and  by  6.  that  I{N2)=  I(Nj)  • 


Excision  of  Relative  Attractors  and  Repellers . 


Definition 

Let  S  be  a  compact  invariant  set  and  let  £  be  a  subset  of  S  . 

Then  S  is  called  an  isolated  invariant  set  relative  to  S  if  there  is  a  compact 

A  A  A  A 

relative  neighborhood  N  of  S  in  S  such  that  S  =  I  ( N )  .  If  this  is 
the  case,  let  A+  *  A+(S\  S )  (resp.  A~  =  A~(£J ,  S)  )  be  the  set  of  points  on 

solutions  in  S  Vs  that  tend  to  £3  in  forward  time  (resp.  backward  time). 


A 

Lemma :  Suppose  N  is  compact,  let  S=I(N),  and  let  S  be  isolated 

4.  4.  S\  —  .  A 

relative  to  S  •  Suppose  either  of  A=A(S,S)  or  A=A(S,S)  is  empty . 

A 

Then  for  all  small  enough  open  neighborhoods,  U  ,  of  S  , 

I  (N\u)  fl  9(N\u)  =  I  (N)  n  8N\A+  UA'US  . 


Proof.  Suppose  A~  is  empty.  Let  N  be  an  isolating  neighborhood  of 

A  A  /\  _ 

S  in  S  .  Since  S  is  the  maximal  Invariant  set  in  N  and  A  is 

A 

non-empty  every  point  in  the  boundary  of  S  n  N  (relative  to  S)  must  leave 
N  in  backward  time .  It  follows  (using  [1],  Ch.  II,  5.10)  that  S  can  be 

A 

written  as  a  disjoint  union  S  U  T  U  S  where  S  is  a  compact  invariant 

~  A 

subset  of  S  and  solutions  in  T  tend  to  S  in  backward  time  and  S  in 


forward  time. 
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Now  let  U  be  any  neighborhood  of  S  with  closure  disjoint  from  S  . 

^  — * 

Then  I(N\u)=  S  and  the  result  follows.  The  argument  when  A  is  empty 

is  the  same  after  a  time  reversal. 


E.  Non-Removable  Critical  Points. 

Lemma :  The  isolated  rest  points  in  I(N^)  are  also  isolated  as  invariant 

sets  relative  to  S  ( Nj)  . 

Let  pa  (Uq,0)  be  a  rest  point  and  let  M  be  the  Jacobean  matrix  of 
(ulfl*u2f2>  at  u0  * 

Assume  det  M  /  0  (then  (uQ,0)  is  isolated). 

Suppose  that  A+  =  A+  (p,  I  ( Nj) )  and  A-  s  A-  ( p,  I  ( Nj) )  are  both 
non-empty  (so  p  cannot  be  "removed"  from  N^  in  the  way  described  in  D.  ). 


Then  : 

1.  If  either  coordinate  of  u0  is  positive,  the  corresponding  diagonal 
entry  of  M  is  negative. 

2 .  If  both  are  positive  then  det  M  >  0  . 

3.  If  both  coordinates  of  uQ  are  zero  then  at  least  one  of  the  diagonal 
entries  is  negative.  Furthermore,  if  the  first  diagonal  entry  is 
positive,  then  the  wave  velocity  9  is  negative,  while  if  the  second 
diagonal  entry  is  positive,  9  is  positive. 

Proof.  As  seen  in  B.  ,  g(u)  =  u^  -  Uj  is  strictly  decreasing  on  non-constant 
solutions  in  I  ( N^)  .  It  follows  (in  general)  that  all  solutions  have  their 
a  -  and  w -limit  sets  in  different  components  of  the  rest  point  set  and  from 
there  that  isolated  rest  points  are  also  isolated  invariant  sets  relative  to  I(N^)  . 


i 
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Now  suppose  p  =  (uQ,  0)  is  a  rest  point  and  let 


(a 

1  • 

1  / 

©  = 

0  ^ 

and  M  = 

f  ° 

i\ 

\c 

dJ 

1°  Hi 

0, 

Then  M  is  the  matrix  of  the  linearized  equation  at  p  and 
det  M  =  det  M  which  has  been  assumed  to  be  non-zero.  It  follows  that 
(u0,0)  is  isolated. 

Suppose  that  X  is  an  eigenvalue  of  M  with  eigenvector 
(4j.42.Tij.Ti2)  •  Let  £  =  (6j,42)  and  n  =  (tij,ti2);  then  ri  =  X|  and 
-M£  -  0ti  =  Xti  .  These  equations  are  equivalent  to: 


(X2  +  0jX  +  a)  |j  +  b|2  =  0 

(12) 

(x2  +  e2x  +  d)42  +  c4j  =  o 


If  A+  is  non-empty  there  is  a  solution  in  I(Nj)  tending  to  uQ  in 
forward  time.  It  must  come  in  along  some  eigenspace  of  M  corresponding  to 
an  eigenvalue  with  non- positive  real  part.  But  it  cannot  oscillate  either,  so 
the  eigenvalue  must  be  real  and  therefore  negative  (since  it  is  not  zero). 

Furthermore,  rij  and  cannot  have  the  same  sign  since  v^  2  0 

and  v2  s  0  on  solutions  in  S(Nj)  .  Since  £  =  Xn  ,  £j  and  £2 
cannot  have  the  same  sign  either  so  it  can  be  assumed  that  4j  a  0  and 

h  s  0  • 

Similarly,  if  k~  is  non-empty,  M  admits  a  real  positive  eigenvalue 
whose  corresponding  eigenvector  satisfies  4j  a  0  and  s  0  . 
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Now  suppose  the  first  coordinate  of  Uq  is  positive.  Then 
b=  8/3u2  ( Uj  fj)  |  p  must  be  negative.  From  the  first  equation  in  (12)  it 
follows  that  for  any  eigenvector,  cannot  be  zero  and  so  the  eigenvalue 

must  have  the  form  A  =  -0^/2  ±  J  ef/4  -(a  +  b^/fij)  • 

Suppose  a  +  b^/tj  is  Positive.  Then  the  real  part  of  any  eigenvalue 
has  the  same  sign  as  -01  =  -6^  or  -0  since  61  >  0  .  Thus  it  is  not 
possible  that  a  +  is  positive  for  both  of  the  eigenvectors.  But 

b£2/£j  Is  non- negative;  it  follows  that  a  is  negative. 

Similarly  if  the  second  coordinate  of  uQ  is  positive  then  d  is 
negative . 

Furthermore,  if  both  coordinates  are  positive,  there  must  be  an 
eigenvector  such  that  both  of  a  +  and  d  +  c^/^  are  negative 

(namely  the  one  corresponding  to  the  eigenvalue  with  sign  different  from  -0  ). 
But  this  easily  gives  ad  -  be  =  det  M  >  0  . 

Finally,  suppose  both  coordinates  of  uQ  are  zero. 

Then  the  solution  in  A+  must  lie  in  U~  while  that  in  A-  must  lie 
in  Uj+  .  The  eigenvalue  for  the  first  must  be  \  =  -9^/2  -  /4  "  d  ana 

that  for  the  second  must  be  -9^/2  ±  n/q^/4  -  a  .  Since  one  of  these  cannot 
have  the  sign  of  -9  ,  one  of  a  or  d  must  be  negative.  Also,  if  d  >  0 
then  the  negative  eigenvalue  (for  the  solution  in  U~  )  has  the  sign  of  -9  , 
so  -9  must  be  negative .  Similarly,  if  a  >  0,  -0  must  be  positive. 

This  concludes  the  proof. 


§  5.  Connection  Triples  For  The  Main  Problem 


A.  Hypothesis  on  the  zero  sets  of  the  F's  . 


It  will  be  assumed  that  and  F2  (and  therefore  fj  and  ) 
have  zero  sets  as  pictured  in  Figure  6. 

There  are,  therefore,  seven 
critical  points ,  of  which  three  are 
attractors  for  the  reaction  equations  , 

(i.e.  in  Cj  ),  namely  (u^,0), 

(0,u2)  and  (0,0).  It  is  assumed 
that  the  Jacobean  matrix,  M  ,  is  non¬ 
degenerate  at  each  critical  point. 

The  squares  are  projections  of  the  se1 
described  in  4C  . 

Let  N'  be  that  about  (0,u2)  and  N”  that  about  (Uj,0). 

The  union  of  the  discs  is  the  projection  in  u -space  of  a  neighborhood, 
U  of  the  remaining  points  which  are  to  be  deleted  from  N2  ( as  in  4E  . ) 

An  isolating  neighborhood,  N3  ,  is  determined  by  N2\U  . 

Now  modify  the  6's  so  that  6^  =  62  ;  of  course  N'  ,  N"  and 
N3  remain  isolating  neighborhoods  throughout. 


B.  A  second  hypothesis  on  the  F's  . 

However,  N3  is  not  the  desired  neighborhood— in  particular,  in 
addition  to  (Uj,0)  and  (0,u2),  it  contains  the  critical  point  (0,0)  . 

It  is  easy  to  see  that  for  some  (unique)  value  of  ,  there  is  a  solution 


of  the  equation 
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(13) 


v 


1 


Vj  =  -SjVj  -  UjfjtUj,  0) 


which  connects  (0,0)  to  (Uj,0)  .  (For  example,  this  could  be  done  by 

continuing  this  system  to  that  in  §2  ).  Furthermore,  the  sign  of  -9. 

,U1 

is  determined  by  that  of  J  f,  (  s ,  0)  ds  ;  namely  ,  it  has  the  same  sign 

0  1 

as  this  integral. 

Similarly,  there  is  a  ( unique)  connection  from  (O,!^)  to  (0,0)  for 
the  equation 


(14) 


v2  =  ‘02v2  "  ^<°'u2)  ' 


^2 

and  it  has  the  same  sign  as  /  f,  (  0  ,  s)  ds  . 

0 

Impose  the  hypothesis  now  that  both  of  the  above  mentioned  integrals 
are  positive.  Then  0^  >  0  >  9j  ,  so  for  any  fixed  values  of  the  6j  (with 
0^  >  0  )  the  point  ( 0 , 0 , 0 , 0)  is  either  an  attractor  or  a  repeller  ( or  both) 
relative  to  I  ( N^)  . 

Using  the  Lemma  of  4.D,  choose  and  delete  from  an  appropriate 

neighborhood  V  of  (0,0, 0,0)  to  obtain  the  last  isolating  neighborhood 
N=N3\V  of  the  triple  (N',N",N). 


C.  Continuation . 

It  is  clear  that  the  F's  can  now  be  modified,  maintaining  the  hypothesis 
in  A  and  B  ,  to  have  the  aspect  in  Figure  7  wherein  three  of  the  critical 
points  have  collapsed.  In  this  situation  it  is  again  clear  that  any  connection 
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from  (0,u,)  to  (0,0)  must  have 

03k 

02  <  0  and  any  one  from  (0,0)  to  \ 

(Uj,0)  must  have  9^  >  0  .  So  again  /r\ 

an  isolating  neighborhood  N  is  obtained  // 

on  deleting  a  neighborhood  of  the  origin.  L— i 

(N'  ,  and  N’ '  are  always  determined  by 

_  Figure  7 

the  squares  about  (u2,0)  and  (0 ,  Uj,) )  . 

Now  it  is  clear  that  the  continuation  can  be  carried  through  the  degenerate 
situation  in  Figure  7  to  the  example  of  §3  with  isolating  neighborhoods  N*  , 
N"  and  N  as  shown  in  Figure  8. 

Now  it  must  be  shown  that  the  - - 

neighborhoods  N'  ,  N”  and  N  LJk. 

determine  the  sets  S' ,  S"  and  S 

of  the  connection  triple  of  the  problem  \  \  . 

® ^  ffl  )■ 

in  §  3.  It  has  already  been  seen  that 

Figure  8 

the  sets  N'  and  N"  determine  the 


Figure  8 


critical  Doints  (0 ,  ,  0 , 0)  and  ( u^,  0 , 0 , 0)  so  it  only  remains  to  see  that 

N  determines  the  correct  set. 

In  order  to  do  this  it  is  sufficient  to  show  that  the  maximal  invariant  set 
in  the  region  indicated  in  Figure  7  is  actually  in  the  invariant  manifold  w^  =  1  . 
To  do  this  consider  the  reaction  vector  field  on  the  two  solid  lines  in  Figure  9 
with  slope  -1  .  In  the  w  coordinates  l 


of  §  3 ,  the  w^  ( =  +  u2 )  component 

of  this  field  is 

4wj  -  3w^  -  Wj3  +  3w2(  1-w3)  .  In 
particular,  it  is  positive  if  Wj  <  1 
and  negative  if  Wj  >  1  .  Thus  the  vector 


\\u 


Figure  9 


SSfi 
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points  into  the  region  between  the  segments  so  long  as  this  region  contains  the 
(dashed)  line  w^  =  1  .  Since  the  negative  of  this  reaction  field  acts  as  the 
"acceleration"  in  equations  (6)  ,  it  is  clear  that  solutions  must  "bounce  off" 
the  lines  as  shown  in  Fig.  8.  (At  this  stage  9j  =  ©2  is  used  of  course.) 

On  collapsing  the  segments  to  the  dashed  line,  it  follows  that  the  maximal 
invariant  set  in  any  such  region  is  contained  in  the  invariant  manifold.  Thus 
the  continuation  of  the  general  problem  to  that  of  §  3  is  completed. 


D.  Connection  triples. 

The  last  step  of  the  proof  of  the  theorem  stated  in  the  introduction  is  to 
show  that  the  triple  ( N' ,  N"  ,  N )  is  always  a  connection  triple  provided  9 
is  large  enough.  This  will  be  true  provided  that  for  each  fixed  6^,  62  (  >  0 ) 

values  0Q  <  0  and  >  0  can  be  chosen  so  that  the  equations  (4)  of  the 
introduction  do  not  admit  connections  from  ( 0 ,  u2 , 0 , 0)  and  ( Uj,  0 , 0 , 0) 
throughout  the  continuation. 

The  idea  here  is  that  if  1 0 }  is  large  then  bounded  solutions  of  (4) 
with  0  ^  u^  ^  K  have  projections  in  u -space  that  are  close  to  solutions  of 
(1)  .  To  see  this,  make  the  transformation  w^  =  -©v^  -  u^  and  multiply  the 
right  hand  side  of  the  transformed  equations  by  -0  .  The  resulting  equations  are 
written  more  conveniently  in  the  following  notation :  let  G  =  ( F^ ,  u2  F2 ) 
and  let  dG  be  the  Jacobean  matrix  with  respect  to  G  .  Also  let 
6  *  =  diag  (  6^  \  &7*)  and  e  =  9  2  .  The  equations  are  then: 

u'  =  eG  +  eW 


w'  =  6  w  -  edG  (w+  G) 
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Since  dG  and  G  are  uniformly  bounded  for  0  s  u^  s  K ,  it  is 
clear  from  these  equations  that  w  j  {bounded  solutions  with  0  ^  u^  s  K} 
is  uniformly  bounded  by  C  £  for  some  constant  C  (  cf .  proof  in  4A  ). 

Therefore ,  in  the  limit  e  —  0  ( |  0  j  —  » )  the  u  components  of  solutions 

satisfy  (1)  (with  a  time  scale  change) .  Of  course  this  is  true  uniformly  over 
the  (compact)  parameter  range  needed  to  deform  the  equation  to  the  standard 
form  in  §  3  . 

Now  the  critical  points  (0,u2)  and  ( u^,  0)  are  attractors  so  it  is 
obvious  that  they  cannot  be  connected  by  a  solution.  Therefore,  for  large 
enough  j  e  |  the  corresponding  points  of  (4)  cannot  be  connected  and  it  follows 
that  (N',N",N)  determines  a  connection  triple  provided  ] 0Q |  and  |0j| 
are  large  enough. 

As  a  consequence,  the  theorem  of  the  introduction  is  proved. 

E.  Concluding  Remarks. 

In  view  of  the  general  lemmas  in  §4  other  arrangements  of  the  zero  sets 
of  the  F's  could  be  considered,  but  the  case  treated  seems  sufficient  for  the 
present. 

More  interesting  is  the  situation  where  the  integrals  in  B.  above,  and 
consequently  the  wave  velocities  0^  and  0-,  in  that  section,  have  the 
opposite  signs  from  those  chosen.  In  this  case  the  deformation  to  the  standard 
problem  is  not  possible  since  the  critical  point  at  Uj  =  u2  =  0  is  not  a  relative 
attractor  or  repeller  in  I(N3)  for  some  parameter  value. 

Without  a  proof,  it  seems  clear  that  the  correct  interpretation  is  the 
following.  In  the  sense  of  D.  of  §1 ,  the  attractor  at  Uj  =  u2  =  0  is 
now  dominant  for  both  of  the  equations  (13)  and  (14)  of  §  5  B.  If  one  imagines 
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the  population  to  be  ( initially)  symmetrically  distributed  in  space  so  that  in 
half  the  space  u^  is  at  its  maximal  stable  level  u^  and  u2  is  zero, 
and  in  the  other  half  u2  =  and  Uj  =  0  then  instead  of  either  of  them 
eventually  dominating ,  one  should  expect  that  the  competition  on  the  mutual 
boundary  drives  both  populations  to  such  a  low  level  that  neither  can  survive. 
Thus  the  competition  plus  the  weakness  of  the  high  level  attractors  relative 
to  the  zero  level  ones  leads  to  mutual  extinction. 

In  reference  [4]  more  mathematical  arguments  are  given  bearing  on 
this  aspect  of  the  problem. 
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The  indices  h(S'),  h(S")  and  h(S',S",S)  for  the  example  in 

§2  can  be  determined  pictorially  as  follows  (cf.  also  [2J). 

First  consider  i>’(  9)  for,  say,  9=0;  this  is  the  left  hand 

(hyperbolic)  rest  point  in  Fig.  A.l. 

The  (  solid)  diamond  shaped  ngion 

about  that  point  is  N  as  well  as 

Nj  .  The  heavily  shaded  boundary 

segments  are  N->  .  Then  [N,  /N,  1 

2  1  2  Figure  A.l 

is  homotopy  equivalent  to  the  pointed 

circle,  S1  .  (The  deformation  squeezes  the  diamond  to  an  arc  and  N2  to 
the  end  points.  Pasting  the  end  points  together,  one  has  the  pointed  circle). 
The  index  of  the  right  hand  end  point  is  seen  to  be  21  in  the  same 

way. 

To  find  h  consider  the  region  N(9)  in  Figure  A. 2  : 
the  curved  arcs  in  the  boundary  are 
orbits  of  the  equation  with  9=0. 

When  0  is  negative  orbits  cross 
out  of  the  region  on  these  arcs  ;  when 
0  is  positive  orbits  cross  into  the 

Figure  A. 2 

region  across  them.  Otherwise,  for 

all  0  ,  orbits  cross  the  straight  boundary  arcs  as  indicated.  Those  where 
they  leave  are  heavily  shaded ;  they  are  in  N2  • 

If  0  s  0  ,  the  top  orbit  segment  is  also  in  N2  ;  for  0^0  the 
bottom  one  is.  In  any  case  N?(0)  has  three  components. 
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Suppose  now  that  9^  >  0  . 
appears  as  In  Fig.  A.  3.  If  0j  >  0 
the  picture  is  as  in  Fig.  A. 4.  In 
either  case ,  the  unstable  manifold 
connects  2  components  of  N2  , 
but  the  pairs  differ.  Let 


Then  the  unstable  manifold  from  S' (0Q) 


N2  =  N2  U  {unstable  manifold  of 
S'(0q)  and  S' (9^)}  .  The  way 
the  components  are  connected  is  then 
measured  by  [N/N2]  .  Up  to 
homotopy  N  is  a  ball  and  N2  ^  H 

consists  of  a  set  which  can  be 

deformed  to  a  point.  (This  can  be  "seen"  with  a  little  patience).  Thus 
[N/N-,]  is  homotopy  equivalent  to  the  one-point  space. 

In  fact,  [N/N2]  is  h.  ( The  definition  of  h  could  have  been 
phrased  in  terms  like  the  above— for  "non-pathological"  situations  at  least; 
however,  the  definition  given  in  §2  is  easier  to  work  with). 

If  both  0Q  and  0^  were  negative  (say)  then  Fig.  A.  3  would  be 
the  relevant  one  at  both  ends  of  [©^ ,  9^]  .  In  this  case  (defining  N2  similarly 
to  the  above)  N2  is  homotopic  to  a  point  and  a  ( disjoint)  circle  in  the  boundary 
of  N  .  Then  h=  [N/N2]  is  22  V  2*  as  stated  in  §2  . 

Perhaps  the  most  precise  way  to  compute  these  indices  is  to  continue 
the  appended  equation  of  the  Lemma  in  2.D.  to  an  even  easier  form  (as  was 
done  for  the  main  problem)  the  detail  would  not  be  very  interesting. 


...  ..... 
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